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ABSTRACT 


This  report  contains  a  description  of  the  models  to  be  used  in 
analyzing  the  capabilities  of  gro\ind-based  sensors  in  determining  the 
mass  of  orbiting  bodies,  model  coefficients,  and  the  justification  for 
their  selection.  Relations  are  derived  for  computing  sensitivity  coeffi¬ 
cients  and  their  coupling  to  mass  variance  for  the  case  of  noisy,  biased 
sensors  (monostatic  and  tri-static  radars.  Baker  -  Nunn  cameras),  and 
for  spherical  and  tumbling  cylindrical  satellites. 
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SECTION  I 


INTRODUCTION 


This  document  is  the  second  in  a  series  of  three  theoretical  reports 
prepared  under  Air  Force  Contract  FI9628-67-COO4I •  The  first  of  the 
series,  Technical  Report  Number  1  of  the  referenced  contract^was  concerned 
with  the  mathematical  models  and  relationships  necessary  to  perform  a 
detailed  maximum-likelihood/minimum-variance  error  analysis  of  the  capability 
of  ground-based  sensors  in  determining  the  mass  of  a  satellite^’'.  It  was 
specialized  to  account  for  the  following  restrictions: 

1)  The  sensors  observed  the  satellite  without  error. 

2)  The  satellite  was  a  sphere  of  5-nieter  diameter. 

3)  All  the  physical  characteristics  of  the  body  except  mass 
were  known  without  error. 

4)  All  the  error  in  the  computed  mass  resulted  from  errors 
and  uncertainties  in  the  knowledge  of  the  orbit-perturbing  forces. 

The  present  report  is  a  companion  document,  extending  the  theoretical 

development  to  remove  the  restrictions  of  perfect  sensor  observations  and 
of  perfect  knowledge  of  the  non-mass  body  characteristics.  In  addition, 
the  body  shapes  are  generalized  from  the  sphere  of  the  preceding  report 
to  include  also  one  of  a  pair  of  cylindrical  objects  of  length  10  meters 
and  diameters  2  and  5  meters,  respectively.  The  cylinders  are  free  to 
tumble  in  a  propeller  motion,  with  the  tumble  axis  perpendicular  to  the 
flight  path. 

Since  the  analysis  in  many  places  depends  upon  the  developments 

Superscript  numerals  denote  entries  in  the  References  Section  of  the  present 
report. 


1 


5 

of  the  previous  report  ,  it  is  suggested  that  the  reader  have  that  document 
available  for  reference . 
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SECTION  II 


PERTURBATIONS  TO  WO-BODY  I-DTION 


A.  GENERAL  EQUATIONS  OF  MOTION 
In  general, 

dx/dt  =  f(x,  t)  +  F(x^  ^  ^  ■t')> 


dj|/dt  =  G;(x,  f,,  u,  t). 


v;here 


x(t) 


X' 
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is  the  state  vector  describing  completely  the  orbit  at  any  time 
the  vector  of  tvro-body  accelerations 


f(x,t) 
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and  F  is  the  vector  of  non-two-body  accelerations 
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vAiere 


*^rag  ~  atmospheric  drag  acceleration  vector, 

Y  ,  =  solar  pressure  acceleration  vector, 

-solar  ^ 

—sun,  moon  =  gravitational  acceleration  vector  due  to  the  sun  and 
moon  taken  as  point  masses, 

F  =  gravitational  acceleration  vector  due  to  the  asphericity 
of  the  earth 

The  vector  ^  can  be  as  large  as  a  six-dimensional  state  vector 
which  describes  rigid-bodjr  rotation,  having  the  Euler  angles  of  the 
satellite  as  its  first  three  elements  and  their  time  derivatives  as  its 
re:naining  three.  The  rotational  excitation  vector  G(x,  u,  t)  re¬ 
flects  the  possible  coupling  of  the  orbit  to  the  rotation,  c.g.,  through 
gravity-gradient  torquing,  and  the  coupling  of  the  dynamic  biases,  as, 
for  example,  drag,  to  these  rotations. 

The  vector  ^  appears  in  F  both  because  the  rotations  present  changing 

cross-sectional  areas  in  r ,  and  'r  ,  ,  and  because  of  the  conscrva- 

— drag  —solar 

tive  exchange  of  energy  between  the  translational  motion  of  the  orbit 
and  the  rotational  motion  of  the  satellite  about  its  center  of  mass. 

That  is,  if  the  satellite  tumble-rate  increases  in  a  conservative  field, 
the  kinetic  energj.^  must  come  from  someidiere,  and  that  "someidiere"  is  the 
translational  energgr  of  the  orbit. 

This  conservative  coupling  is  investigated  in  Appendix  I  of  this 
report,  where  it  is  shovmi  that  it  is  indeed  possible  to  de-orbit  a  satel¬ 
lite  in  a  non-dissipative  environment  using  only  torque  devices.  Ilov.'evcr, 
the  minimum  dc-orbiting  time  is  shovrn  to  be  on  the  order  of  10"'’  years, 
so  this  conservative  coupling  of  ^  to  x  can  be,  and  is  in  (2.3)j  ignored. 
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B.  TUSSLING  MODE 


The  mode  of  tu:i±iling  for  non-spherical  satellites  is  restricted 
under  this  contract  to  be  a  propellor-like  notion,  vdth  the  axis  of  ro¬ 
tation  directed  along  the  satellite's  velocity  vector.  Since,  also,  the 
non-spherical  satellites  are  designated  to  be  cylinders,  it  is  clear  that 
the  same  cross-sectional  area  of  a  satellite  is  alira.ys  presented  to  the 
atmospheric  drag  mechanism.  As  far  as  is  concerned,  then,  it  is 

not  a  function  of  the  rotations ^(t) . 

Since  the  sun  and  earth  see  changing  aspect  areas  for  this  mode  of 
tumbling,  the  direct  and  earth-reflected  solar  pressure  vector  iSj 

hovrever,  a  function  of  ^(t). 

T;:o  things  must  now  be  done:  the  precise  d^manics  of  IgCt)  must  be 
established  or  assumed,  and  the  functional  dependence  of  *r^^^  upon 
^(t)  must  be  derived. 

1.  Dotation  D^-namics 

In  order  tliat  the  tumbling  plane  turn  vdth  the  satellite  velocity’’ 

A 

vector,  i.e. ,  in  order  that  the  velocity  unit  vector  i^  alv;ays  be  normal 
to  the  tumble  plane,  as  shovm  in  Figure  1,  belov^,  torques  internal  to 
the  satellite  must  be  ass'omed.  For  slnplicity,  it  vdll  be  assu;;ied  that 

these  interval  torques  are  such  that  the  tiuible  rate,  say  is  a  con- 

A 

stant  relative  to  the  normal  i  to  the  orbital  plane.  Because  of  the  i.on- 

n  ^ 

A 

tvro-body  effects,  i^  vdll  move  in  inertial  space,  and  so  the  tumble  rate 
vdll  not  be  inertially  constant.  Since  the  nutation  rate  of  the  orbital 
plane  is  relatively  small,  the  satellite  torques  necessary  to  explain 
this  inertially  accelerating  motion  vrould  also  be  small.  In  any  event, 
the  i.iotion  assumed  is  as  valid  as  anj'  other  once  the  propeller  motion  is 
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FIGURE  1.  TUMBLING  GEOMETRY 
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presxiraedj  and  the  results  of  the  error  analysis  ;vill  at  least  be  indi¬ 
cative  of  the  importance  of  accurate  signature  data. 

The  tumble  model  is^  in  algebraic  terms, 

t(t)  + ,  (2.4) 

’idth  the  rotation  state  ?  defined  as  the  2-vector 

?  [1]  . 

« 

If  the  errors  in  the  constants  and  vrere  to  be  modeled  (a  covariance 
;.;atri:c  for  them  introduced  in  the  estimation  equations),  then  G  in  (2.2) 
vrould  have  to  be  explicitly"  vn'itten  to  express  the  dynamics  of  constaiit- 

velocity.  bliile  this  is  very'  simple  to  do,  ire  iri.ll  talie  the  errors  in 

# 

^  ^  and  to  be  unmodeled  errors  and  find,  byr  the  quasi-experii'.iental 
means  discussed  in  Section  IV. D  of  Keference  5,  the  sensitivity  coeffici¬ 
ents  This  iri.ll  indicate  the  relative  importance  of  signature  in- 

fonnation  iri.thout  depending  on  the  reasonableness  of  any  assumed  value 
for  a  covariance  matrix  on 
2.  Solar  Pressure  as  a  Function  of  Tumble 

The  details  of  how  *f'  ,  functionally  depends  on  S  are  derived  in 

Appendix  II  of  this  report.  In  summary',  the  results  are  as  follows. 
Define; 

=  specular  reflectivity  of  satellite, 

=  diffuse  reflectivity  of  satellite, 
q  =  Earth  albedo, 

X  =  ratio  of  earth  radius  to  distance  between  satellite 
and  earth  center, 

=  arc  cos  X  , 

I  =  solar  constant, 

c  =  speed  of  light. 
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i  =  unit  vector  from  center  of  earth  to  center  of  satellite, 
r 

1.  =  unit  velocity  vector  of  satellite 
r 

^  =  unit  vector  from  center  of  sun  to  earth, 

ss 

A  A 

cos  a  =  -(i  .  i  ), 

r  S3 


|a|  <rr/2  -  % 
tt/2  -  3^  4  |c  I  ^  rr/2  + 
rr/2  +  |c  I  rr 


1(1  +  >qa2) 
(1-1- 


|o|  4  n/2  -  3^ 

t\/2  -  \  ^\c.\  ^  rr/2  + 

n/2  +  f:  jc  I  <:  n 


=(-.0417  +  .5431><)/3, 

a^  ={:0444  -  3.17(4  -.77)^  +  .0045  (7-. 77)  sin  [u.3(>'-.77)n]}/3 
%s  ""  ®2  ^  “  se®"^^  -  e“'^y  (2  -  sy)]  /2 

\  ^{^ss  +  (a/2)  C  s  +  1  -  s(l  +  87)‘^J  }  cos  a 

f  CClA-sin  a)^  +  (4-sin  a)^]  -  (1-4^)-^^ ""  ~i  sin  a 

^  (l+>i^-  2>'  sin  cy^^-  ^  ^ 


t”  =  -4  +  9.34 
y  =  (c-n/2)/B^ 

f  1  j  0, 

(-1,  7>0, 

d  -  3.7  +  59(>-.77)^  , 

A  A 

Sin©  =  i  •  i.  , 
r _ r 

cos  0  =  \/l-sir?  , 

--  (i^X  ij.)/  cos©  ,  the  normal  to  the  orbital  plane. 
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H-> 


n 


A  A  A  A 

i  =  i. X  i  ,  the  vector  normal  to  i  that  defines  the  tumble  plane 
p  r  n 

(see  Figure  1), 

_  =  sin  ?  1  +  cosZ  i 

3  ^  P  n 

A  A 

COS  a  =  i  •  i„  , 
s  ss  3  ' 


COS  a  =1 
r  r 


"3  ^ 


cos  0  ---  (p^  cos  +  p^  cos  +  p. 


Y' 


A  =  ttR  ,  the  cross-sectional  area  of  a  sphere  of  radius  R  , 
s  s 

2 

-  nR^  ,  the  area  of  the  base  of  a  cylinder  of  radius 
A^  -  2LR  ,  cross-sectional  area  normal  to  the  base  of  a  cylinder  of 

JL  N./ 

length  L. 


Then  the  solar  radiation  acceleration  is,  for  a 


Spherical  Satellite: 


Cylindrical  Satellite: 


y  ^  =  (I  /c)  (A,/m) 

-solar  s'  '  T' 


(2.5) 


1/'-  I  cos  0  I 

+  TTk^/6  +  (1  +  k^3)  I  sin  0  1  ]  [p^  +  p^  igg] 

+  [2  (VV  ^^r  1+  -  (V3)k^  I  sin  0  I 

-  ]  [  Pi  +  P2  “s  ]  S  ]  (2.6) 

Taking  the  "ballistic"  coefficient  with  respect  to  solar  pressure 


to  be 


u^= 


(Inojj/c)  (A/ra)  ,  spherical  case, 
(Ij^on/°^  (Ai/m)  ,  cylindrical  case. 


(2.7) 
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where 


I  =  I  (R  /r  )^  , 

s  non  es'  es 


vriiere  R  is  the  mean  earth-sun  distance  and  r  is  the  actual  earth-sun 
es  es 

distance,  it  is  clear  from  (2.5)  and  (2.6)  that  the  ^  vector  that  repre¬ 
sents  simply  the  terns  in  (2.5)  or  (2.6),  respectively, 

enclosed  in  the  large  square  brackets  times  (^g^/  ^es^  ' 


y  =  (R  /r  )‘ 
-  es'  es 


terns  in  large  square  brackets  in  (2.5)  or  (2 


.6)1 


(2.S) 
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SECTION  III 


ESTIMATION  AND  PREDICTION 


A.  PRELIMINARY  DETAILS 

Consider  equations  (4.1)  and  (4.2)  from  Reference  5,  suitably 
modified  to  account  for  the  tumbling  of  a  non-spherical  satellite: 

d(Ax)/dt  =  f(x,  t)  -  f(x°,  t)  +  F(x,  ,  u,  t),  (3.1) 

Az  =  h(x,  ,  b,  t)  -  h(x°,  0,  0,  t) ,  (3.2) 

where  Ax(t)  is  the  Encke  variation  from  the  two-body  orbit  x°(t) 
to  the  actual  orbit  x(t)>  T  is  the  vector  of  two-body  accelerations, 
and  F  is  the  vector  of  non-two-body  accelerations. 


F(x,|,u,  t)  = 


't, 

-drag 


+ 


^olar 


+ 


^  — 

0 

p - - 

+ 

0 

p — 

sviiy  moon 

r 

—0 

developed  in  detail  in  that  reference  and  in  Section  II  of  this  report. 

The  vector  u  contained  in  F  represents  the  dynamic  biases  affecting 
the  satellite.  The  first  three  of  its  components, 

v^(t)  = 

are  "ballistic"  coefficients  which  will  be  actively  estimated.  The 
remaining  three  of  its  components  account  for  uncertainties  in  the 
oblateness  forces  F^  due  to  errors  in  x  and  in  the  geopotential  coef¬ 
ficients  J.  These  are  defined  by  the  3-vector 

=4(x,  J)  -  F^(x,  J), 
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where  x  ~  best  estimates  of  x(t)  and  J  at  time 


t  given  all  data  upto  and  including  time  To  first-order  accuracy, 


(3.3) 


where  r  and  Ar  are  the  first  three  components  (the  position  components)  of 


X  and  Ax,  respectively.  The  hats  over  the  partial  derivatives  indicate  that 
they  are  evaluated  around  r(t  and  u(t/t^_^). 

The  2- vector  ^  (t)  contained  in  F  describes  the  rotations  of  the  cylin¬ 
drical  satellites  to  be  considered  in  this  phase  of  the  study.  Since  there 
are  no  very  obvious  covariance  matrices  applicable  to  (t),  it  will  be 
treated  as  an  unmodeled  parameter.  In  keeping  with  the  discussion  on  sensi¬ 
tivity  coefficients  in  Section  IV. D  of  Reference  5,  the  study  will  there¬ 
fore  seek  the  sensitivity  matrix  in  order  to  ascertain  the  impor¬ 

tance  of  errors  in  our  knowledge  of  tumble  position  and  rate.  A  detailed 
analysis  of  the  tumbling  mode  called  out  in  this  study  is  performed  in 
Section  II  of  the  present  report. 


In  (3.2), 


z  =  h(x,  ,  b  ,  t) 

defines  the  observations  ^(t)  on  the  orbit,  where 


(3.4) 


=  zero-mean  Gaussian  white  noise 


b  =  constant  but  unknown  biases 


A^(t)  “  ^  ~  the  difference  between  the  actual  noisy  obser¬ 


vation  and  the  noise-free  observation  of  the  two-body 


part  of  the  motion. 
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Now,  we  can  define  the  error  vector 


e  (t/t.  )  =  a  (t)  -  a  (t/t,  )  for  any  vector  a  (t).  In  the  special 
^  ^  cases  where  t  assvines  discrete  values  t^,  the  notation 

will  be  simplified  to  (n/k). 


V/ith  the  ^-vector  y(t)  defined  as  the  column  of  variables  which  are  to 
be  actively  estimated,  i.e.. 


z(t)  = 


(3.5) 


the  difference  equation  relating  the  one-step  prediction  errors  to  the 
prior  in-step  estimation  errors  is  (equation  (4.31) >  Reference  5) 

(3.6) 

,  ey(n/n-l)  =  $y(n,  n-1)  j^(n-l/n-l)  +  —  (n,  n-1)  e^  (n-l/n-1)  +  ^(n-t), 

“2  2 

where 


$y.(n,  n-1)  = 


$(n,  n-1) 


d  X 

8j<n,  n-1) 


/ 


T 

n 


-  t 


n-1’ 


$ (n,n-l) 


=  the  two-body  transition  matrix  for  the  states  x  (see 


Appendix  IV,  Reference  5)j 
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3(a),  y (a) 


dr^rag  /^^]_  ^^olar  evaluated  around  the 

best  state  estimate  x(t/s), 


dv 


—  (n,  n-1)  =  r"  |(t 


s)  ds 


0  0 

0  0 
0  0 
0  0 
0 
1 

0  0  0 
0  0  0 
0 


o 
0 
0 
1 

0  1 
0  0 


0 


0 


I. 

a(n-l)  =y  s) 


n-1 


Wj^(s) 

0 

0 


Wj(t)  =  zero-mean  Gaussian  white  noise  with  power  (2/'^^) 

per  unit  double  bandwidth  in  rad. /sec. 

Finally  the  following  definitions  for  covariance  matrices  will  be  needed. 

For  any  vectors  a(t)  and  b(t), 

Cov  [a(t),  b(t^  =  E[^a(t)  ^t)]  -  E[a(t^  E[^b'^(t'^  , 

where  E(  )  denotes  mathematical  expectation  and  (  denotes  matrix 
transpose^ 

and 

P^t.(nA)  =  Ccv[a(n)  -  ^(nA),  ^n)  -  ^(nA^  ,  [eJnA).  e^,(nA)] 

The  covariances  on  the  en^dronment  are 


Qyy  (n)  =  Cov  u(n), 


(n)  =  Cov  Vn), 


\b 


=  Cov  b 
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with 


and  all  cross-variances  among  w,  and  b  identically  zero. 

B.  MBUMOM- VARIANCE  DERIVATION 
1.  Fundamental  Expressions 

The  derivations  depend  upon  three  basic  expressions:  the  one-step 
prediction  error  equation  already  presented  as  (3*6),  a  linearized  obser¬ 
vation  equation  developed  from  (3*2),  and  the  in-step  estimation  equation 
which  develops  from  the  linear  regression  solution  to  the  minimum-variance 
problem. 

The  linear  regression  solution  is  characterized  as 


y(n/n)  =  y(n/n-l)  +  By.(n)  [  Az(n)  -  Az(n/n-l)] 
where  B  (n)  is  a  gain  matrix  to  be  determined  and  where,  from  (3-2), 

y 


(3.7) 


Az(n)  -  A^(n/n-l)  =  hfx,  b,tj|  -  ©7^,  t^ 


3.8 


After  a  first-order  expansion  around  x(n/n-l),  the  latter  can  be  written 
as 


Az(n)  -  Az(n/n-l)  =  Hy.(n)^  ^(n)  -  l(n/n-l)]]  +  H»|^(n)ij(n) 

+  Hj^(n)b 


3.9 
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where  H  (n)  = 


I 

H(n)  I 


0  0  0 

•  •  • 

•  •  • 

•  •  • 


10  0  0 


H(n)  = 


ah 

dx 


(n)  =  ah 


H^(n)  =  ah 
db 


X  =  x(n/n-l) 


X  =  x(n/n-l) 


X  =  x(n/n-l) 


Subtracting  ^n)  from  both  sides  of  (3.7)  and  employing  (3-9), 
we  find 


ey.(n/n)  = 


fl  -  B  (n)H  (n)1 
L  y'  7  J 


^(n/n-1)  -  By(n)H^(n)  q  (n) 


-  By(n)H^(n)b. 


(3.10) 
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2.  Coiaputation  of  Optimal  Gain  B  (n) 


As  discussed  in  Reference  5>  the  minimum- variance  gain  matrix  is 
found  by  setting 


trace  |a  (n/n)  a|  =  min. 


(3.11) 


No  matter  idiat  the  form  of  the  iveighting  matrix  A  (in  Reference  5  vre 
happened  to  assume  it  diagonal),  the  minimum  estimation  variance  obtains 
from  the  solution  to 


dP  (n/n)  =  0. 

yy  ^ 

Since,  from  (3.10jl 


P  (n/n)  =  [l-E  H  ]  P  (n/n-l) 
jy  '  7  7-^  yy 

+  B  H,  R,  ,  R  ^B  -  [l-B  H  ]  P  , 
y  b  t)b  1)  y  y  y  yb 


[l-B  K  ]  '^  +  B  H  R„„(n)  H»,'^B 

y  y-*  y  ^  ^  y 

(n/n-l)H, B  ^  -  B  H,  P  (n/n-l)  [l-B  H  ] 

W  /  b  y  y  y  /  /  L  y  yJ 


T 


(3.12) 


vdiere  B  and  the  H's  have  argument  (n),  it  can  be  shown  in  the  manner  of 

y 

Reference  5  that  the  optimal  gain  is 


By(n)  =  [Pyy(n/n-l)Hy'^(n)  +  Py^( n/n-l )Hj^'^(n)]  [Hy(n)Pyy( n/n-l )Hy'^(n) 
+  (n)R,^,^  (n)H^'^(n)  + 

+  Hy(n)Py.j^(  n/n-l  )H^'^(n)  +  Hj^(n)Py,^(  n/n-l  )Hy.'^(n)^ 


(3.13) 
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Note  the  covariauice  laatrijc  P  ,  (n/r-1)  that  appears  here.  If,  as  in  the 

j  ^ 

perfect  sensor  case,  we.  are  to  generate  the  estimation  covariances  recur¬ 
sively,  two  nev/  expressions  must  be  derived  in  addition  to  those  that  ;;ere 
needed  for  perfect  sensors;  P^^^(n/n-l)  as  a  function  of  P^^Cn-l/n-l) ,  arid 
P^^(n-l/n-l)  as  a  function  of  B^^(n-l)  and  Py.^(n-l/n-2). 

3.  Piediction  Covariance  Computations 


Directly  from  the  perfect  sensor  case  , 


1) 


+  -DP  (n-l/n-1) 


(n.  -DP^3(n-Vn-l)J 

2  V 


C^(n-l)  ^  (3,13) 


r  A  ^ 

dF  I 

’yj(n/n-l)  =  yn,  n-1)  +  —  (n,  n-1)  ^  (n-l)J  P^jCn-l/n-l) 


A 

dr  dF 

+  ^  (n,  n-1)  ^  (n-1) 


(3.U) 


;\(here  P  ^  =  Cov  (aJ),  the  covariance  of  the  geopotential  coefficient 

<J  c; 

errors, 

T 

Multiplying  (3-6)  by  b  and  taking  the  expected  value,  ;ve  compute 
the  additional  prediction  covariance  needed  for  this  real  sensor  case: 
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(3.15) 


P  ,  (n/n-l)  =  8  (n,n-l)  +  (n,n-l)  (n-1)  P  ,  (n-l/n-1). 

L  /  av^  a^r  -I 


4.  In-Step  Covariance  Computations 

Upon  substituting  the  optimal  gain  B  (n)  defined  by  (3.13)  into 
equation  (3.12),  xife  can  rearrange  terms  and  make  the  appropriate  can¬ 
cellations  to  obtain 

P^,(n/n)  =  [l  -  B^^(n)  H^,(n)]  P^/ n/n-l)  -  Bjn)  H^(n)  Py^^( n/n-l) 


yy 


y 


yy 


(3.16) 


This  is  to  be  compared  to  (4.41)  in  Fteference  5. 


lie  can  post-multiply  (3. 10)  by  e  (n/n)  =  y  (n),  use  the  definition 
for  y^(n)  provided  by  (3.3)^  ^^nd  take  the  ecxpectation: 


Pj^^(n/n)  =  [l  -  By(n)  H^Cn)  ]  [p^( n/n-l) 

]- 


^  j. 

— o  (rj  + 

ay 


A  ip  A  p 

P  -(n/n-l)  o  (n)  -  B  (n)  H,  (n)  P  (n/n-l)  o  (n) 


(3.17) 


xdiere  the  argument  (n)  on  the  partial  derivatives  indicates  evaluation 


at  the  point  x(n/n),  and 

ap 

— c 

ay 


-o  = 


ap  I  coo  I  000 
■T^  1 000  I  000 
-  1 000  I  000 

L.  I  I  . 


Equation  (3.1?)  is  to  be  compared  with  (4.42)  of  Reference  5. 


The  matrix  P  (n/n)  does  not  differ  from  its  perfect-sensor  form: 
^2^2 
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A  ^  T 

^  (n/n)  =  ^  (n)P^^(n/n)  ^  (n) 
2  2 


A  A  fp 

+  ^  M  +]l(n/n)  +  II^(n/n), 


dJ 


where 


A  A 

dF  d? 

Tl(n/n)  =  (n)  Pj.j(n/n)  ^  (n). 


(: 


However,  a  nev/  matrix  need  be  computed  for  the  right  side  of  (3.15). 

T 

Post-multiplying  (3. 10)  by  b  ,  we  can  calculate 

(3.19) 

to  complete  the  error-analysis  recursions. 

The  filtering  process  is  the  same  as  in  the  perfect-sensor  case, 
except  that  the  P^^matrix  recursions  are  new  and  the  microscopic  de¬ 
tails  of  the  other  recursion  relationships  reflect  the  observation- 
noise  and  bias  statistics. 


.18) 
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SECTION  IV 


SENSOR  MODELS 


Three  basic  sensor  types  are  to  be  studied  under  the  subject  con¬ 
tract  for  their  performance  in  mass  determination.  These  are  (l)  mono¬ 
static  radars ,  vrLth  observation  vector 

Range 

_  Range-Rate 
— '  '  Azimuth 

Elevation 

(2)  tristatic  radars .  vri.th  observation  vector 

=  j  b,  t), 

and  (3)j  Baker-Nunn  cameras ,  vdth  observation  vector 

fj.\  Right  Ascension  _  ,  /  v  ^  \ 

[Declination  J  t). 

The  detailed  functional  relationship  between  the  observations  and 
the  states  are  presented  in  Appendix  III,  where  the  partial  derivatives 
with  respect  to  states,  noise,  and  biases  are  developed.  Note  that  all 
the  sensors  are  affected  by  time  biases,  and  the  Baker-Nunn  cameras  are 
affected  by  a-c  time  fluctuations  as  well. 


Range 

Range  Difference 
Range  Difference 
Range-Rate 

Range-Rate  Difference 
Range-Rate  Difference 
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APPENDK  I 


THE  COUPLING  OF  SATELLITE  TUMBLING  TO  ITS 
ORBITAL  MOTION  IN  A  CONSERVATIVE  FIELD 

1  2  S 

Various  writers  *  have  considered  the  problem  of  the  motion  of  a 
non-symmetrical  satellite  in  the  gravitational  field  of  a  major  planet. 

Although  the  satellite  may  nominally  be  travelling  in  a  stable  orbit,  there 
is  obviously  an  interaction  between  the  equations  vdiich  govern  the  rotation?! 
motion  of  the  satellite  about  its  center  of  mass  and  the  equations  of  orbital 
motion.  In  some  cases,  interest  in  this  problem  centers  about  the  advantageous 
use  of  this  interaction  for  purposes  of  gravity  gradient  stabilization;  in 
other  instances,  there  is  concern  about  the  possibly  detrimental  effects  of 
this  coupling  such  as  excessive  tumbling  of  the  satellite,  or  even  significant 
distortion  of  the  orbital  path. 

It  is  one  of  the  purposes  of  this  appendix  to  show  that,  in  almost  any 
practical  situation,  these  detrimental  effects  are  entirely  insignificant. 
Furthermore,  the  perturbation  of  the  orbit  due  to  the  lack  of  spherical 
symmetry  of  the  satellite  can  almost  always  be  ignored,  and  the  rotational 
(tiuriblin^  equations  can  be  solved  with  insignificant  error  under  the  assumption 
that  the  center  of  mass  of  the  satellite  is  moving  in  a  pre-computed  orbit.  Using 
this  assumption,  the  tumbling  equations  for  the  satellite  will  be  derived  in 
a  form  which  is  useful  for  computational  purposes. 

In  the  following  derivation,  it  will  be  assumed  that  the  satellite  is 
a  body  of  revolution,  or  at  least  that  two  of  the  three  principal  moments  of 
inertia  are  equal.  The  results,  however,  leave  little  doubt  that  the  above 
claims  are  equally  valid  for  a  body  having  all  three  principal  moments  of 


inertia  unequal. 


In  Reference  1,  it  is  shovm  that  the  gravitational  torque  N  exerted  on 
a  body  in  an  "inverse  square  law"  gravity  potential  field  is  given  by 

v/here  (x)  denotes  the  vector  cross  product,  and 


gravitational  constant; 
r  --  distance  betv/een  centers; 

K  A  A  A 

i^  =  unit  radial  vector  expressed  in  an  i^,  i^,  i,,  coordinate 
system  v/hich  is  fixed  in  the  body; 

AAA 

J  =  inertia  matrix  of  the  body  idth  respect  to  the  i^,  i^,  i, 
coordinates, 


A  A  A 

If  the  i^,  i^,  i^  coordinates  are  directed  along  the  principal  inertial  axes 
of  thie  body,  then  J  is  a  diagonal  matrix  having  diagonal  elements  equal  to 
the  principal  moments  of  inertia. 

Let  w  denote  the  angular  velocity  vector  v\rLth  components  equal  to  the 
angular  velocities  about  the  principal  axes,  Euler's  equations  of  motion 
describing  the  rotation  of  the  body  have,  in  this  coordinate  system,  the 
well-knovai  form  ^ 


Jw  r  (0  X  J(o  =  N, 


where  the  dot  denotes  differentiation  with  respect  to  time.  In  the  case  of  a 
body  of  revolution,  we  take 


J  = 


I 

0 

0 


0 

0 

I-AI 


(I.l) 


(1.2) 


G 

I 

0 
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(  L3) 


where  AI  is  positive  for  long  cigar-shaped  bodies,  negative  for  saucer-shaped 
bodies. 

Let  the  orthogonal  matrix  S  represent  a  rotation  through  an  angle  p 
about  the  symmetric  principal  axis  of  the  body,  and  set  w  =  SX,  or 


cosp 

-sinp 

0 


sinp 

cosp 

0 


0 

0 

1 


X. 


If  the  angle  p  is  allowed  to  vary  with  time,  thai 


•  •  * 

w  =  SX  +  SX,  vAiere 


s  = 


sinp 

cosp 

0 


-cosp 

sinp 

0 


0 

0  . 

0 


We  also  introduce  the  skew  33rmmetric  matrlxjl,  defined  as 


Jl  = 


0 


so  that  we  can  write  w  x  Jw  =nJco.  With  this  notation,  (1.2)  becomes 


J(SX  +  SX)  +  ilJSX  =  N. 

T  * 

Since  J  obviously  coimirutes  with  S,  S  and  S,  this  latter  equation  can  be 
put  into  the  form 

JX  +  UX  =  s'^N, 

T  *  T/s. 

with  L  =  SS  +  SliS.  If  X^,  X^,  X^  are  the  components  of  the  vector  X, 


( 1.4) 
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we  find 


0 


L  = 


-(p-^3) 


0 


K 


1 


and 


LJ\  =  (ip  -  AI\^) 


(1.5) 


T  ^  ^ 

With  respect  to  equation  (I.l),  note  that  S  (i^  x  is  simply  the 

A  A 

vector  i^  X  Ji^  expressed  in  the  rotated  coordinate  system  having  unit  vectors 

A  TA 

=  S  (n  =  1,  2,  3).  this  system,  the  inertia  matrix  J  still  has 


the  form  (1,3);  hence, 


=  Ii  -  Al(i  •  j„;j„ 


in  ^3^  coordinates,  and 


J3)(ir  ^  V 


= 


-3 


d 

r 

-(t 


Jl)(^r*  V 


(1.6) 
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Equations  (l.4)}  (l.5)  and  (1.6)  show  that  =  0,  or  X^Ct)  =  'K^(O) . 
If  we  now  let  =  (i^*  n  =  1,  2,  3j  we  can  write  the  equations  of 
rotational  motion  in  the  following  simple  form. 


d 

dt 


^2"’3 
“^1^3 J 


Let  us  now  relate  the  equation  (1.7)  to  a  coordinate  system, 

(n,p,q),  having  unit  vectors  with  ineH-ially  fixed  directions,  and  an 
origin  moving  with  the  satellite.  Assume  that  the  satellite  is  travel¬ 
ling  in  a  simple  elliptic  orbit.  (This  will  certainly  be  valid  for 
sufficiently  small  time  intervals;  moreover,  we  intend  to  show  that 
the  perturbation  of  the  orbit  due  to  tumbling  is  actually  negligible  over 
a  period  of  many  orbits.)  Define  the  orthogonal  unit  vectors  and 

as  follows: 


ip  =  1^  at  the  instant  the  satellite  is  passing  through  perigee; 


i  =  unit  vector  in  the  direction  of  instantaneous  velocity  at  perigee; 


i^  =  unit  vector  normal  to  the  orbital  plane,  in  the  direction  of  the  orbit 
momentum  vector. 


If  V  is  the  true  anomaly,  then  in  this  fixed-direction  system  the 
vector  ij.  has  the  representation 


A 

1  = 
r 


cos  V 
sin  V 
0 


(1.7) 


(1.8) 
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p 


FIGURE  2.  TUl'IBLING  GEOI-ETRY 
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Figure  2  shows  the  orientation  of  the  satellite  in  the  fixed  system,  in 
terms  of  both  the  fixed  body  coordinates  and  the  rotating 

^1^  ^2^  ^3  *ngle  p  has  been  chosen  in  such  a  way  that  the 

unit  vector  is  in  the  (q,n)  plane.  It  is  clear  that 


a 

X  =  Y  sin  a 

P  +  Y  cos  a 

In  (n,p,q)  coordinates,  we  have 


^1 


angular  velocities  along  . 


0 

sin  a 

cos  a 

A 

A 

cos  Y 

^  Jg  - 

-cos  a  sin  Y 

,  ^3  = 

sin  a  sin  Y 

sin  Y 

cos  a  cos  Y  _ 

-sin  a  cos  Y_ 

(l.9) 


(LIO) 


and  we  find 


r^  =  cos  Y  sin  v 


A  A 

^r  *  ^1 


r^  =  sin  a  cos  v  -  cos  a  sin  Y  sin  v 


r  '^2 


A  A 

r^  =  cos  a  cos  v  +  sin  a  sin  Y  sin  v  =  i^  • 


(  lOl) 


If  a  and  e  are,  respectively,  the  semi-major  axis  and  the  eccentricity  of 
the  orbit,  it  is  well  known  that 


1  +  e  cos  V  * 


(1.12) 
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Noting  finally  that  (AlX--Ip)=  (Al-I)X.  +  lY  cos  a,  it  follows  from  equations 

3  3 

(l*9)j  (I. 11)  and  (1.12)  that  the  tumbling  equations  (1.7)  can  be  written 
in  terms  of  the  physically  meaningful  angles  a,  Y  and  v,  and  the  constant 
It  is  of  particular  importance  that,  in  their  final  form,  these  equations  are 
entirely  independent  of  the  angle  p,  which  in  the  case  of  a  body  of  revolution 
is  an  unobservable  quantity. 

The  kinetic  energy  of  rotation  is  given  by 


T  =  (l/a)^”^!^  =  (l/2)(o'^(SJs’^)(o  =  (l/2)(o'^(JSs'^)(o  =  (l/2)u'^J(o.  (1.13) 


Differentiation  with  respect  to  T  produces 

•  T  • 

T  =  X  JX, 


(1.14) 


which  may  be  regarded  as  either  the  rate  of  increase  of  the  rotational  energy 

of  the  satellite,  or  the  rate  of  decrease  of  energy  in  the  orbit.  Premultiplying 
T  T 

(1.4)  by  X  ,  and  noting  that  X  LJX  =  0,  we  find 


•  T  T 
T  =  X  S  N 


=  4M  ''I  ^3 

where  the  latter  expression  follows  from  (l*6). 

In  order  to  justify  the  previous  statements  concerning  the  tiimbling  rates 
and  the  perturbation  of  the  orbit  due  to  tumbling,  it  would  be  sufficient  to 
show  that  the  change  in  T  over  a  period  of  several  orbits  is  negligible  both 
with  respect  to  the  total  orbital  energy  and  the  implied  change  in  tumbling 
rate.  The  direct  computation  of  the  energy  or  the  energy  rate,  however,  would 


(1. 15) 
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essentially  require  the  solution  of  the  nonlinear  tumbling  equations  (I. 7)>  and 
this  problem  is  extremely  difficult  if  not  intractable. 

Consider,  instead,  the  follovd.ng  hypothetical  situation.  Suppose  that 
the  satellite  is  equipped  with  a  set  of  torque  rockets  which  can  produce 
arbitrary  rotations,  but  no  translational  motion,  of  the  satellite.  Suppose, 
in  addition,  that  these  rockets  are  either  pre-programmed  or  controlled  by  an 
intelligent  being  in  such  a  way  that  the  total  decrease  in  orbital  energy 
over  a  pre-specified  time  interval  is  maximized.  It  is  assumed  that  the  rockets 
are  capable  of  producing  both  continuous  and  impulsive  outputs,  so  that  there 
are  no  inertial  limitations  on  the  possible  rotational  motions  of  the  satellite. 

Since  the  equations  (1.7)  no  longer  govern  the  satellite  rotations,  (1.15) 
does  not  represent  the  rate  of  increase  of  rotational  energy,  but  (1.15)  is 
still  the  correct  expression  for  the  rate  of  decrease  of  orbital  energy.  Further¬ 
more,  the  maximum  decrease  in  orbital  energy  under  the  artificial  motion 
described  above  is  an  upper  bound  for  both  the  actual  decrease  in  orbital 
energy  and  the  actual  increase  in  rotational  energy  under  any  natural  motion 
of  the  satellite. 

Since  we  are  still  assuming  motion  in  a  nominal  elliptic  orbit,  we  can 
treat  T,  a  andy  as  functions  of  timet  or,  equivalently,  as  functions  of  the 
true  anomaly  v.  Noting  equations  (1.9)  and  (1.12),  we  rewrite  (1. 15)  as 

dT  _  3j^(l+e  cos  v)^  AI  r 
dv  a^(l-e^)^ 

where  r^,  and  r^  are  still  given  by  (I.ll).  Our  plan  is  now  to  integrate 
(A. 16)  between  specified  values  of  v,  and  to  maximize  the  resulting  expression 


[r^  ^  “  ^1  —1 

dv  dv  J 


(1.16) 
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with  respect  to  the  unknown  functions  a(v)  and  y(v).  For  simplicity,  we  will 
work  with  the  absolute  value  |aI].  Without  displaying  all  of  the  cumbersome 
details,  we  state  that  the  function Y(v)  which  produces  a  maximum  for  this 
integral  is  simply  "'((v)  =  tt/2.  This  can  be  obtained  either  as  a  solution  of 
the  Euler-Lagrange  equation^  which  is  associated  with  this  variational  problem, 
or  by  recognizing  that  x  Ji^)  is  largest  when  the  angular  velocity 

about  an  axis  in  the  direction  of  x  ^  maximum,  and  this  occurs 

vdienY  =  n/2. 

Setting  y (v)  =  n/2  in  (1.16),  and  substituting  for  r^  and  r^,  we  obtain 

dT  _  3/4  ( 1+e  cos  v)^  |aI|  sin  2  (a-v)  ^ 

dv  2a^(l-e^)^  dv 


Integrating  this  expression  between  limits  v^  and  v^,  we  are  left  with  a  simple 
calculus  of  variations  problem  in  the  single  unknown  function  a(v).  The 
associated  Euler-Lagrange  equation  is  algebraic,  and  leads  to  the  result 


tan2(a-v) 


-2(l+e  cos  v) 

3e  sin  V 


(1.18) 


From  (1.18),  we  compute 


da(v)  _  I  ,  _ 3e(&bcos  v) _ 

p  P  , 

dv  (3e  sin  v)  +  4(l+e  cos  v) 


and 


sin2(o.-v)  = 


2(1  He  cos  v) 


1^(3  e  sin  v)^  f4(l+e  cos  v)^  ] 
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Clearly  vrtien  e  <k/l>  we  niust  choose  the  positive  sign  for  sin2(a-v),  since 
the  expression  for  da/dv  is  positive  for  all  v  in  this  case.  ActuaJLly,  a 
closer  analysis  shows  that  the  positive  sign  should  be  chosen  for  all  e<l, 
but  the  sinqjler  observation  is  sufficient  for  practical  purposes. 

At  the  endpoints  of  the  interval  of  integration,  if  the  angle  a(v)  does 
not  satisfy  equation  (l.l8),  an  impulse  from  our  hypothetical  rockets  can 
be  used  to  bring  the  satellite  into  position.  Inspection  of  (I.l?)  leads 
to  the  conclusion  that,  if  the  integration  is  performed  over  an  integer 
number  of  orbits,  the  maximum  contribution  to  the  integral  due  to  endpoint 
effects  is  given  by 


'  a3(l-e)3 


(I 


Substituting  the  expressions  for  sin2(a-v)  and  da/dv  into  (A.I7)  and 
integrating  over  a  single  orbit,  we  find 


AT  <  3>^n|Al|  F(e), 
orbit  a^(l-e^)^ 

where 

N  _  1  r  2(1+8  QOS  v)^[C3e  ain  v  )^  +  4(l+e  cos  v)^  +  3e(e+  cos  v)]  dv 
^  Jq  [(3e  sin  v)^  +  k  (1+e  cos 

The  following  table  shows  the  results  of  a  numerical  integration  for  the 

2 

function  F(e),  as  well  as  a  tabulation  of  the  values  of  1  +  2e  . 


19) 
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e 

F(e) 

1+2  e^ 

0 

1.000 

1.000 

0.1 

1.021 

1.020 

0.3 

1.184 

1.180 

0.5 

1.503 

1.500 

0.7 

1.970 

1.980 

1.0 

2.939 

3.000 

Hence,  we  have  the  very  simple  estimate 


AT  <  3iU  n|Al|  a+2e^) 

orbit  ~  3  r-i  2x3 

ar  (1-e 


(I. -20) 


Combining  the  results  (1.19)  and  (l.20),  we  obtain  the  following  upper 
bound  for  the  decrease  in  orbital  energy,  and  the  corresponding  increase  in 
tumbling  energy,  over  a  period  of  n  orbits. 


|AT(n)l<4^' 


(l+e)~'  +  nn(l+2e^) 


(1.21) 


Since  the  total  energy  of  the  orbit  is,  initially,  h  =  m/(/2a,  vdiere 
m  is  the  mass  of  the  satellite,  the  fractional  change  in  orbital  energy  per 
orbit  (neglecting  endpoint  effects)  is  no  greater  than 


^  ^  IaiI  (l+2e^) 

h  a”^  m  (1-e'^)'' 


(1.22) 
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For  the  extreme  case  of  a  cylindrical  satellite  having  length  L  and  radius 
R«L,  this  becomes 


^  <  n  (h\  (l+2e^) 

h  ~  2  la/  (l-e^)2 


(1.23) 


n 

For  an  earth  orbit,  we  must  have  a(l-e)>  2  x  lO'ft.  Also,  for  0<  e<0.8, 

2  3 

the  inequality  (l+2e  )/(l+e)^(l-e) < 2  is  valid.  Hence,  over  the  stated  range 
for  e,  we  have  Ah/h<10~^^^,  where  L  is  the  length,  in  feet,  of  the  cylindrical 
satellite. 

Assume  now  that  the  decrease  in  orbital  energy  over  a  single  orbit 
corresponds  to  an  increase  of  I(Ao))  /2  units  of  energy  due  to  a  change  in  the 
angular  velocity  of  the  cylindrical  satellite  about  some  non-syrametric 
principal  axis.  In  this  case,  we  must  have 


Mi  LLjgs!!  .  (1.24, 

[a(l-e)]  2  I  (n-e)3 

Noting  that  |aI|/I<1  for  any  cylinder  and  (l+2e^)/(l+e)^^  1  for  e>0,  we 
16  3  2 

take  1.4  X  10  ft.  "^/sec.  and  the  above  inequality  for  a(l-e)  to  obtain 

(Aw)^<33  X  10~^  rad.^/sec.^  (1*25) 

for  a  complete  orbit.  This  estimate  is  extremely  conservative,  since  the 
satellite  cannot  tumble  end  over  end  and  simultaneously  maintain  a  configiiration 
which  causes  the  orbital  energy  to  decrease  at  each  instant.  The  main  value 
of  the  estimate  (1.25)  is  in  demonstrating  that,  over  the  small  portion  of  an 
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orbit  which  lies  within  the  range  of  a  single  tracking  station,  the  net 
change  in  angular  velocity  due  to  gravitational  effects  is  rather  small. 
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APPENDIX  II 


PERTURBATIONS  ON  CYLINDRICAL  SATELLITES 
DUE  TO  SOLAR  RADIATION  PRESSURE 


A.  GENERAL 

There  are  three  components  of  solar  radiation  force  acting  on  a 
satellite: 

(1)  That  due  to  energy  impinging  directly  from  the  sun  and  re¬ 
flected  from  the  earth. 

(2)  That  due  to  the  reaction  of  specularlj’"  reflecting  radiation 
from  the  satellite. 

(3)  That  due  to  the  reaction  of  diffusely  scattering  radiation 
from  the  satellite. 

As  noted  in  Appendix  I  of  Reference  5?  if  the  incident  radiation 
approaches  from  the  direction  shox^m  in  Figure  3  below,  then  the 
pressure  exerted  normal  to  is 

E  =  (Ig/c)  =  p  (Il.l) 

v/here  I  is  the  irradiance  and  c  is  the  speed  of  light.  V/hen  £  impinges 
s 

upon  a  differential  surface  area  dA,  the  differential  force  due  to  the 
incident  flux  is 

dF.  ='1.  (p  •  d/)  ='1.  (I  /c)  (1.  .  dA),  (II. 2) 

— xn  1  -  1  s'  1  — 


vz-hile  the  differential  force  due  to  the  flux  being  reflected  back  from 
the  surface  is 


dF 

— refl 


A 

=  -i 


refl 


(I  /c)  (i 


refl 


dA), 


(II.3) 
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FIGURE  3  .  TUMBLING  COI'T^IGURATIOM 
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A 

where  i 


is  the  unit  vector  in  the  direction  of  the  reflected  flux  and 


refl 

I  is  the  irradiance  in  that  direction, 
r 

The  net  solar  force  due  to  direct  solar  radiation  is  the  integral 
of  (II. 2)  plus  (II. 3). 

Now,  the  total  pressure  £  at  the  satellite  consists  of  a  part  di¬ 
rected  from  the  sun  and  a  part  reflected  from  the  earth.  For  simpli- 

A 

city,  this  pressure  vector  xvill  be  resolved  into  components  in  the  i^ 

A 

and  i^  directions,  i.e.,  radial  from  earth  and  sun-to- satellite,  re- 
G  5 

spectively.  This  pressure  vector  can  be  written  as 


£=  (I/c)  (Pi  t 

5 


ivhere,  to  a  good  approximation. 


(II. 4) 


Pi  = 


P2  = 


^  4q  (a^  +  a^)  cos  a 
Xq  a 

r 

0 


f 


(1  +  Xqa^) 
(1  +  Xqa^^) 


given  that 


1  a  I  <  rr/2  - 

n/2  -  4  |c  I  -^r  tt/2  + 

rr/2  +  Bj^  ^  |c  I  tt 


|a|  ^  rr/2  -  B. 


rr/2  -  B^  4|a|  ^  rr/2  +  B^ 
n/2  +  ^  |a  I  <  rr 


q  =  Earth  albedo, 

X  =  ratio  of  earth  radius  to  distance  betxveen  satellite 

and  center  of  earth 
=  arc  cos  X, 
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»(-.0417  +  .543lX)/3, 

*|.0444  -  3.17(X  -.77)^  +  .0045  (7-. 77)  sin  [l4.3(^-. 77)n]jr/3 

^ss  ""  ^2  ®  -  e~'^j  (2  +  sy)  ]  /2 

a  =la  +  (a,/2)  [  s  +  1  -  s(l  +  sy)*^]]  cos  a 
r  I  ss  i.  ^ 

+  f  CClA-sir.  +  (X-sin  a)^]  -  (1-X^)^^^  T  sin 

^  (1+X^-  2>-  sin  X  j  6 


^  =  -4  +  9.3X 

y  =  (a-TT/2)/^^ 


d=  3.7  +  59(X-.77)^  . 

B.  IJET  Bi'DIATION  FORCES 
1.  Incident  Radiation 

IJeglecting  the  reflection  effects,  the  impinging  energy  i.'ill 

A 

produce  a  net  force  in  the  i^  direction  defined  by  the  integration  of 
(II. 2)  over  the  surface  of  the  satellite: 


— 

F.  = 
—in 

s 

dF.  =  r 
—in  ^ 

[  i.  .  dA  + 

- 1 

•H 

<*H 

(II. 6) 

■'.d'lere 

9 

--  ,  the  area  of  the  base  of  a  cylinder  of  radius  I.  , 

A,  =  2Ln  ,  cross-sectional  area  normal  to  the  base  of  a  cylinder 
1  c^ 

length  L. 

A 

Taking  the  unit  vector  i„  to  be  normal  to  the  base  of  the  cylinder,  (see 

Figure  3),  and  defining 

A  A 

COS  0  =  -(i  •  i^). 
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where  ^  "*■ 

;/e  find  that  (II.  6)  beconies- 


^  I  ^  I  +  I  sin  0  |  ). 


(II. 7) 


2.  Specularly  Reflected  Radiation 

Consider  a  raj'-  of  light  -which  is  reflected  from  an  increment 
dA  of  surface  area  according  to  Snell's  La-w.  If  the  ray  impinges  at 
an  angle  0  i  ^  ^  -  Tr/2  i^rith  the  normal,  then  it  ^^all  leave  at  an  angle 
^  ^  -idth  the  normal,  as  sho-ivn  in  Figure  4  belcjw.  If  the  specular 


FIGURE  4.  SFECUL/Ji  RLFLECTIOII 


reflection  coefficient  is  defined  to  be  k  ,  then  the  reflection  ir- 
- r 

radiance  is  I^  =  I^,  and  the  integration  of  (II. 3)  jm.elds 


-T^efl 


=  -k  p 
r^ 


A 

i 


refl  ^^refl 


dA)  + 


(  * 

side 


refl  ^^refl 


Ro'w,  over  the  base  of  the  cj'-linder. 


dA  =  i^  sign  (cos  0)  dA  , 


dA) 

(II. S) 
(II. 9) 


v/here 


sign  (x)- 


F  • ::: 


'+1  ,  x>0 

the  sign  function  being  necessary  to  distinguish  -i\rhether  the  top  base 
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A 

of  the  bottom  base  of  the  cylinder  is  exposed  to  the  light.  If  i^ 
is  the  normal  to  dA,  it  is  clear  from  Figure  4  that 

Wl^  ^i  2(cos  (II.IO) 

But  from  Figure  3  we  see  that 


cos 


y'n  ""  I  cos  0  I  , 


so  that  (II. 10)  becomes 


r  /\ 

irefi  =  +  2  |cos0|  sign(cos  ^  )  i^ 


=  +  2cos  P  • 


Hence,  over  the  base  area, 


(TI.ll) 


i  ■  dA  =  (  -  I  cos  2 1  cos  0  I  )  dA 

reil  —  I  I 

=  I  cos  <b\  dA. 


(11.12) 


Combining  (II. 11)  and  ( 11.12)  in  the  first  integral  in  (II.8),  we  have 

/  (Irell  '  “A)  =  \  +  2  %)• 

(11.13) 

Now,  integrating  over  the  side  surface  of  the  cylinder,  we  note  that  we 
can  decompose  p  into  two  parts:  One  part  tangential  to  the  longitudinal 
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A 

‘3 


FIGURE  5 .  SPECULAR  REFIECTION  FROM  CYLINDER  WALL 
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axis  and  one  part  normal  to  that  axis.  Figure  5  provides  the  details. 
Only  the  normal  part  of  that  having  direction  cosine  -  sin  i^, 
appears  in  the  dot  product 

A  A  A 


A 

i 


refl 


"a  =  -"i  •  ^A’ 


so  that 


A 

i 


A 


refl 


i.  =  sin  0  sin  a,  0  ■£  0  <  tt , 


(II.U) 


A  A 

where  a  is  the  polar  reference  angle  to  i^.  The  full  unit  vector 
also  is  needed.  Referral  to  Figure  6,  below,  where  the  origin  has  been 
moved  to  the  surface  increment  d/i  of  the  previous  figure,  provides 

i  .  =  (sin  0)(sin  2a)  i,  -  (sin  0)(cos  2a)  i^ 


refl 


-  (cos  0)  i„ 


(11.15) 


Substituting  ( 11.14)  and  ( 11.15)  into  the  second  integral  in  (II.8), 
we  find 


/ 


A  ,A  . 

i  -T  (i  „  •  dA)  = 

refl  '  refl  - 


side 


A)=  /"[(si 


sin  0)(sin  2a) 


-  (sin  0)(cos  2a)  i^  -  (cos  0)  i^  J  (sin  0)(sin  a)  R^Lda, 
=  -  A^  (-  ^  i^  +  (cos  0)(sin  0)  i^), 

=  -  - 3 -  |_4(cos  0)  i,  +  i^J  , 


(11.16) 


0  <  0  ^  n. 
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/s 

*3 


FIGURE  6.  GEOMETRY  FOR  DIRECTION 

A 

COSINES  OF  i  ,  SPECULAR  CASE, 
refl’ 
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The  net  specular  reflection  force  now  follows  from  putting  (II. 13) 


and  (11.16)  into  (II.8): 

-refl  "  "  P  “  (A^/3)Isin  0|1  i^ 

+  2  cos  0  ^  A^j  cos  0|  -  (2A^/3)lsin 

(11.17) 

3.  Diffusely  Reflected  Radiation 

For  the  case  of  diffuse  radiation,  the  reflected  energy  obeys  Lambert's 
Law.  As  developed  in  Appendix  I  of  Reference  5,  the  differential  reaction 
force  is  normal  to  any  increment  of  area  dA: 

'^-refl  ~  ^A  ^^d^s/c^  ^^i  *  * 

(11.18) 

where  is  the  di ffuse  reflection  coefficient.  Again 


— refl 


For  the  base, 


/  / 

'^K  ''  A  .  , 

b  side 


dF  „ 
—refl 


=  -  (2/3)(k^Ig/c)  A^  cos  0  1^. 


(11.19) 
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For  the  side,  (11.14)  can  be  used  to  obtain 


/. 


dF 

-refl 


side 


=  -  (2/3)(k^Ig/c)(sin  0)  J'  j(cos  a)i 
+  (sin  “)i2j 


=  “A^(n/6)(k^Ig/c)(sin  0)i2, 


=  A^(TT/6)(k^yc)(i.  +  cos  0 


Hence,  for  diffuse  reflection, 

4efl  "  ^^d  [irrl^/2)\  +{(nA^/2)  -  2A^)(cos  0)1^]  .  (11.20) 

C.  FINAL  RESULTS 

Combining  (II.?),  (Il.l?),  and  (II. 20),  and  upon  dividing  by  the 
mass  m,  we  obtain  finally,  after  using  (II.4)  to  substitute  for  p 
(or  1^), 

^solar  =  ^  [  [( W  ^  • 

+  nk^/6  +  (1  +  kj3)  I  sin  0  |  ]  [p^  +  p^  J 

[?  ~  I  I 

-  ^d/6  ]  [  Pi  ^  +  P2  .  (IT  21) 
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where 


cos  a 


s 


A  A 

COS  a  =  i  * 

r  r  ^ 


y 


y 


cos  0  =  cos  cos  a^)/ 


The  last  identity  follows  from  the  fact  that 


cos  0  =  _i  •  i 

i 

and  from  (II.4),  which  yields 


A  A  /y  _ _ _ _ 

il  =  (  Pl  Ij,  +  P2  ijs  )//p^2  +  ^2 


US 


APFENDIX  III 


OBSERVATION  EQUATIONS 

The  notation  used  in  this  appendix  in  some  cases  contradicts  the 
definitions  for  the  same  symbols  used  elsewhere  in  this  report  and  in 
Reference  5 •  This  occurs  because  particular  symbols  are  commonly  ac¬ 
cepted  to  mean  one  thing  in  orbital  mechanics ,  and  something  entirely 
different  in  radar  technology. 


The  symbol  definitions  within  this  appendix  which  are  to  be  inter¬ 
preted  locally  only  are  the  following: 

R  =  range 
R  =  range-rate, 

A  =  Azimuth, 

E  =  elevation, 
a  =  right  ascension 
S  =  declination, 

r^=  station  location  in  inertial  geocentric  coordinates, 
range  differences,  i  =  1,2, 

* 

AR^=  range-rate  differences, 

^  ,0*  =  defined  forms  of  station  latitude  and  longitude. 


A.  MOKOSTATIC  RADARS 


The  monostatic  radars  to  be  evaluated  in  this  study  measure  ranro 
R,  range-rate  R,  azimuth  f ,  and  elevation  E.  Given  the  state  vector  ^(t) 
and  the  station  location  vector 

r 

— 3 

in  inertial  coordinates,  and  given  the  inertial  spin  rate  of  the  earth 
'jQ  can  \n:ito  explicitly  the  entries  of  the  nonostatic  ol^servation 
vector 


z(t)  - 


I'' 

J_j 


h  (x,  2[,  b,  t). 


(-1.1) 


These  are 

H  =  1/  (x  -  Xg)^  +  (y  -  y^)^  +  (z  -  , 


R  = 


A  =  tan 


(1/r)  [  (x  -  Xg)(x  +  w^y^)  +  (y  -  y^Xy  -  +  (2 

-1 


-X 


(III. 2) 
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where 


x' '  ' 

sin  cos 

o' 

sin  <p  sin  o' 

-  cos  cp 

1 

X 

1 

ca^  ^ 

1  t  1 

y 

-  sin  o' 

cos  o' 

0 

1 

CO 

_z'" 

cos  (5  cos 

o' 

cos  (p  sin  o' 

sin<p 

1 

_ 1 

sin¥>  =z^  //  +  +  -  ^ 


sin  9'  -7^  // 


(latitude) 

(longitude) 


Then 


dR/dx 

dR/dx 


(1/R) 

(1/R) 


,  z  -  Zg,  0,  0,  oj  , 

R(:!^  y  )  -  R(x-x^),  R(y-<j  x^)  -  R(y-y„) 


X  -  X  ,  y  -  y 
•'s 


e"  s 


e  3 


R 


R 


Rz  -  R(z-z  ) 

_ __s 

R 


,  x-x  ,  y-y  ,  z-z 

s'  *^3'  3 


dk/dx  = 


:  (2''')^  ^(-x'''  sin  0 '  -  y ' ' '  cos  o'  sin  <p  ) ,  (x ' ' '  cos  0  ' 


dE/3x 


V  r2  -  (z"')2  LI 


-  y ' ' '  sin  0  '  sin  ^  ) ,  y ' ' '  cos  ^  ,  0 ,  0 ,  0  J  , 

j”  |cos  o'  cos  (p  -  (z  ' '  '/R^)(x-Xg)^  , 

|sin  o'  cos  {5  -  (z" '/R^)(y-yg)J  , 

|sin  ^  -  (z" '/R^)(z-Zg)|  ,  0,  0,  oj  , 

(III  .3) 


51 


\irtiich  sire  to  be  used  in 


H  = 


dR/dx 
aP./dx 
dA/ dx 
dE/ dx 


(III. 4) 


The  noise  matrix  H  is  the  4X4  identity  matrix 


H  = 


3  matrix  H, 
b 

is  the  4 

X  8  array 

1  1 
1  1 

T(t) 

'  dh 

1  ' 

dr 

1  ^ 

1  I 

O 

1 

1 

—  * 

(III. 5) 


(III.(^. 


associates  ;d.th  the  bias  8-vector 


b  = 


Range  Bias 
Range-Rate  Bias 
Aziiiiuth  Bias 
Elevation  Bias 

Station  Location  Bias 
y  Station  ^cation  Bias 
Station  Location  Bias 
Clock  Bias 


(III. 7) 


Nov^,  the  station  location  biases  are  given  in  an  earth-fixed  coor¬ 
dinate  system  and  are  constant  in  that  system.  Defining  the  earth-fixed 
1  111 

system  r  =  column  (x.  ,  y  ,  z  )  used  in  Appendioo  III  of  Reference  ’.;o 
have  that 

r  =  T(t)  r^  =  PNG  rt 
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where  P,  N,  and  G  are  orthogonal  matrices  that  account  respectively  for 
the  earth's  precession,  nutation,  and  rotation  about  its  polar  axis. 


Note  that  —  is  coiriposed  of  the  rows 

dr 

— s 

~  -3R/ 3r, 

3R/3r^ 

II 

1 - 1 

[R(y-WgX)  - 

3A/ 3r^ 

=  -3A/ or. 

SZ/  £r^ 

3 

=  -3E/6r, 

1/here  ^  is  the  ^-ele'.'rient 
dv 

row  vector  in 

3x 

.y)  -  R(x-x  )]  , 


(III. 8) 


r~  *  -1 

l^or  j  3r  J 


The  timing-bias  part  of  (III. 6),  3^dt,  is  a  column  array  of  the 
following  scalar  partials: 


an/dt  - 

(ap/ar)  r 

+ 

(aR/ar^) 

r 

— s 

aR/at  = 

( aft/  ox)  i 

+ 

(aE/ar^) 

r 

— s 

dA/at  = 

(aA/ar)  i 

4- 

(aA/ar^) 

r 

— s 

aE/at  - 

( aE/ ar)  r 

■h 

(aE/ar„) 

r 

— s 

’.-■here 


-s 


0) 

e 


(III. 9) 


53 


B.  TRI-STATIC  RADARS 


The  tri-static  radars  provide  a  range  and  range-rate  measurement 
R  and  ft  from  one  site,  and  the  differences  AR^  =  R-R^  ,  AR^  =  ^“^2  ^ 
AR^  =  H  -  R^  ,  AR^  =  R  -  F.2  calculated  with  respect  to  the  ranges  and 
range-rates  R^,  R^  and  R^  to  two  other  sites.  The  observation 
vector  is  now 


z(t)  = 


R 

AP,^ 

AR. 

R 

AR, 


=  h  (x,  b,  t), 


(III. 10) 


If  the  three  sites  are  located  respectively  at 


s 

\l‘ 

'"s2 

^s 

'  -si 

^sl 

'  -s2 

^"s2 

CO 

_ 1 

-"si. 

_^s2_ 

then  R  and  R  are  the  same  as  in  the  monostatic  case,  whereas 


AR 


=  R  -  \/(x-x  .)^  +  (y-y  .)^  +  (z-z 

L  V  SI  •'si  SI 


AR.  =  R  - 
1 


_AP  L 

i 

+  (3-7gi)  ^]’ 


1  =  1,2. 


( -Til.  11) 
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The  state  partials  dR/dx  and  aR/dx  are  clearly  the  same  as  before,  but 


d(AR.  )/3x  =  1 

^  R(R-AR^) 


jj(x-Xg)  (R-AR^)  -  (x-Xg^)  r]  , 

[(y-yg)  ' 

[(z-Sg)  (R-AR^)  -  (z-Zg^)  -]  ^ 

0,  0,  0  J  , 


a(af:.)/dx  ^ 


ftCx-x^)  _  +  (R-AR.)  (x-x^.)  j 


R-AR. 

1 


(R-AR. )' 
^  1 


L  „2  R-AR.  x2  J 


R" 


(R-AR. )' 


j-Ri  -  fi.(z-Zg) 


R 


-X-X  x-x  . 

IS-  SI 


R-AR. 

1 


(R-AR^)  (z-Zg^) 
(R-AR^)2 


1  .  rlz^  -  ]  3  r  ^  ^si-i 

J  L  R  R-AR.  J  L  R  R-AR.  J 


.  R 

i  =  1,2. 


R-AR. 

1 


y-; 

R 


z-z 


z-z  . 

SI 

1. 

1 


(III. 12) 


These  arc  rov:  entries  in 


H  = 


aR/  ax 
a(AP,^)/ax 
a(AP2)/ax 
dT./dx 
a(Aft^)/ ax 
a(AR2)/ax 
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The  noise  niatrix  is  now  a  6  x  6  identitv  matrix,  and  the  bia: 


ratrix  I!^  is  the  6  x  16  matrix 


I 

'  T(t)  1 

ah  1 

T(t)  \ 

ah 

1  1 

1 

aAr  .  1 

-s2  ' 

^  - 

v.'hero 


ah 

ar 

— s 


aR/ ar^ 

a(AE^)/arg 

d{^R^)/dr^ 

ap/ar„ 

aCAR^y/ar^ 


oh 


a(AP^)/aArg^ 

0 

0 

a(AR^)/aAr^^ 

0 


ah 


C 

0 


a(AE2V  ^^Is2 


a(AR2)/aAr 


^2 


(III. 14) 
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dh 

3t 


dF/et 

d(AR^)/dt 

a(AR2)/dt 

d^dt 

d(  'R^)/dt 

a(AR2)/at 


Iloi;,  the  surveys  for  the  three  stations  ars.  such  that  the  uastcr 
station  at  r^^  is  located  by  geodetic  survey.  Its  position  has  survey 

The  slave  stations  arc 

urveyed  relative  to  the  raster: 


biases  in  the  earth- fisced  system  (X"^,  Z^). 


r  . 
— 


r 

— s 


'  r  . 

—SI 


i  1,  2, 


(III. 15) 


nd.th  thw  relative  biases  on  the  relative  nosition  v.  cters  Ar  .  inch.- 

—51 

pvcnJont  of  and  much  s^.'^aller  than  the  absolute  survey  bias  or.  th^  vast  .r 
station. 


The  bias  vector  is  the  16  x  1  array 


r  R 

AR^ 
AR. 


.  1 
^7 


Range  Bias 
Range  Bias 

Rai:ige  Bias 

hvaiige-i.ate  Bias 

iiaiige— iLat  e  liias 


haiigc— i-^at  e  n  1.0.3 


X  lister  Ztatioxi  Location  Bias 

y^  Xaster  J 'cation  Location  Bias 

Ixastor  otation  Location  Bias 
Ax^  Slave  1  Lolatlve  Location  Bias 

clave  1  Relative  Location  Bias 
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,  1 


olave  1  r.clativc  Location  Lias 
Slave  2  relative  Location  Bias 
Slave  2  Iiclativo  Location  Lias 
Slave  2  LelatiVv^  Location  Bias 


]lock  Biac 


Upon  putting  (111,15)  into  (111,11),  deflniiig 


..  -  n-en. , 
1  1 


1  1 


and  talcing  derivatives,  \\rc  find,  for  the  range  variables, 


cL/dr^  -5L/3r  , 


;  5 


)/er^  ar./3r,  c(Ar. )/e(^r,S 

5(-rn/3(^r,p  -  (l/i;,)  fcx-x  ,)  ,  (7-7=.n  .  (z-z^nl 

uhcrc  the  second  equation  results  frora  expanding  AR^  into  h-I 
noting  that  -dR./3r^  “  3(^1.  )/'^(hr  .)•  the  range  rates, 

oh/Or^  -  as  given  in  (III.S)  , 
d(Ar.J/£r  -•=  3R/ar.  +  £(A",  )/3(^r^J 

-L.  S  o  1  ^  1 

d(L.2.  )/d(.\r  .)  =  _1  rfh.  (iHo  y)  -  R.  (x-x  .)], 

1  '  '  —SI  - -  I  *-  1  o''  '  l'  SI  -*  » 

Fr.  (:r-ui  x)  -  R.  (y-y  .  )1  , 

Li'-  c 

[R.t  -  R.(z-z  .)  ]  1 

1  si'  -•  J 


R.' 

1 


(III. 16) 


(III. 17) 

L,  and 


(III. 18) 
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The  tiiue  derivatives  for  6^dt  arc 


u:,/5t  =  (an/^r)  r  +  (aR/ar^)  , 
a(iR^)/at  =  [a(AR^)/5r']  t  +  [a(AR^)/ar_]  r^ 

+  [a(AR.)/a(ir^.)]  Ar^.  , 

aR/at  =  (au./a:c)  x  (ar./ar^)  r„  , 

a(Ar.^)/at  =  [a(AR^)/ax  j  x  f[a(Zirt^)/ar^  ]  r^ 

+  [a  ( AR^ )  /a  ( Ar  )  ]  Ar , ^ 


(III. 19) 


C.  BAKER— IIUh’K  G/iZjR^3 


The  Bale cr-I.'uriTi  ca.'acras  record  right 
the  target.  Thesv.  arc,  rcspectivelj'. 


ascension  and  declination  of 


(III. 20) 


•vhere  }f  is  the  right  ascension  of  Greemdch,  the  "Greenvdeh  Hour  Angle." 
The  observation  vector  is 


The 

partial  i 

derivatives 

are  as 

foUous : 

at*  /ax  - 

1 

r  ^-^’s  , 

1 

r\  r\  r\ 

9  Vy 

9 

2 

sec  a  1 

o 

a  <5 /ax  = 

-(“-“s’ 

\  > 

(y-Tg) 

O  o 

,  (z-z  ;  -  R 

S  9 

R-'cos  S  1 

L 

(z-z  ) 

0,  0,  0  , 
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for 


dot /ax 

a  c5/dx 


(III. 22) 


The  noise  rnatriie,  since  a  fluctuating  tine  error  is  assuned  for  this  case^ 
is  nov;  the  2x3  array 


K 


n 


1  0  I  ao<./3t' 
0  1  1  dS/dt_ 


(III. 23) 


;;hcrc 

-  od  /ar^  -  a«./0r  ^ 

—  d  6  / ar^  ~  a  5/ar  , 

and 

aoc/at  =  (aoc/ar)  r  +  (aoi/ar^)  -  ay /at 
a  <5  /at  “  ( a  <5  /ar )  r  +  ( 5  /ar  )  r  , 


(III. 24) 


v.tiere,  to  a  high  degree  of  accuracy  for  orbital  durations  on  the  order 
of  3’’earSj  o  ZT/at  =  o)^  =  const.  These  are  also  used  in  the  2x6  natrix 


H 


b 


1 

0 


0 

1 


cl 


ar 


I(t)  j 

I 


ad /at 
a  <5 /at 


\diich  corresponds  to  the  bias  vector 


night  Ascension  3ias 
Declination  Eias 
x^  Station  Location  Bias 
y^  Station  Location  JBias 
z"^  Station  Location  Bias 
Clock  Bias 


(III. 25) 
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